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We discuss basic features and new developments in recently proposed induced gauge theory jl| solvable in any 
number of dimensions in the limit of infinite number of colours. Its geometrical (string) picture is clarified, using 
planar graph expansion of the corresponding matrix model. New analytical approach is proposed for this theory 
which is based on its equivalence to an effective two-matrix model. It is shown on some particular examples how 
the approach works. This approach may be applicable to a wide class of matrix models with tree-like quadratic 
couplings of matrices. 
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1. Introduction 

The recently proposed D-dimensional induced 
gauge theory solvable in the large N limit ||l| has 
virtually produced some ambitious hopes to pro- 
vide the long searched Master Field (MF) solu- 
tion for QCD. This means that we would find a 
set of field variables that become classical in this 
limit. In principle, we know such a set: Wilson 
loop functionals, which satisfy classical loop equa- 
tions (Makeenko-Migdal equations). But owing 
to their complicated functional character and es- 
sential non-locality they do not advance us very 
much on the way to an analytical solution of mul- 
ticolour QCD. 

Another idea which so far was developed sep- 
arately in field theory, was the so called induced 
gauge theory. Brifly the idea of induction is the 
following: let us start from a Lagrangian of a field 
theory containing only massive gauged matter of 
some kind, with no addition of any kinetic terms 
for the corresponding gauge field. If we integrate 
out the matter we obtain a (generally) non-local 
effective action for only gauge fields. If one then 



increases the mass of matter field one obtains, in 
principle, the l/(mass) expansion of the effective 
action in terms of local operators. 

Of course, things are not as simple as they seem 
at first sight. We have to satisfy a few conditions 
to induce a needed gauge theory: the sign of the 
coupling of an induced interaction should be cor- 
rect, the power-like divergences should be absent 
(to be left with only a logarithmic dependence of 
the induced gauge coupling on the cut-off), non- 
perturbative (e.g. in all orders of l/(mass) ex- 
pansion) locality of the induced action, e.c. 

Surprisingly enough, one can formulate a gauge 
theory of this type: on the one hand, this induces 
(at least naively, in the logarithmic approxima- 
tion) the Yang-Mills gauge coupling, and on the 
other hand, it is solvable by means of the large N 
matrix model technique. The MF appears to be 
local in this case, which immediately raises some 
doubts on the possibility for this theory to serve 
as a right QCD realization in 4 dimensions. And 
indeed, as will be clear from the following, the 
theory obeys symmetry too much (with respect to 
the centre of the gauge group) to describe a real 



2 



QCD §, §. It is doubtful that this symmetry 
can be broken spontaneously, unless we make the 
model more complicated (and thus unsolvable). 

Nevertheless, it seems to be interesting in itself 
that a solvable non-abelian gauge theory of any 
kind exists in physical dimension. 

First of all, the search for an MF for real QCD 
along these lines should not be necessarily a hope- 
less enterprise. Some attempts of this kind are on 
the way 

Secondly, the corresponding D-dimensional 
matrix model might be a meaningful string theory 
with an infinite mass spectrum of physical states. 
In particular, the analytical methods proposed for 
the investigation of this model can be used for 
the solution of any multi-matrix model with the 
tree-like quadratic form of couplings between the 
matrices (which includes, say, Q-component Potts 
models on random graphs) Q. Moreover, the in- 
duced gauge model under consideration is itself 
the model embedded into a tree, rather than into 
a real D-dimensional space, at least in the large-N 
limit 0. For large N the space almost falls out 
from the model and its dimensionality appears 
only through the number of nearest neighbours 
of the D-dimensional lattice. In the next orders 
of 1/N, the space restores little by little, but the 
trivial space structure of the lowest order also sig- 
nals that it is not a QCD realization. 

On the other hand, as a matrix model, this 
theory obeys all the features of some string the- 
ory, which means that it is a theory of an ex- 
tended object from the point of view of the inter- 
nal (transversal) modes of the world sheet given 
by the corresponding planar graphs. Therefore, 
amid the tree nature of the model in the embed- 
ding space, it can produce an infinite and appro- 
priately scaled spectrum of physical states due to 
these transversal modes. So, it is worth trying it. 

An analytic approach of the model was sug- 
gested in 1^, and a considerable progress was 
achieved in the papers Q, where the critical be- 
haviour of the model was established. In the pa- 
pers (see also [|l^,|ll[), the case of Gaussian 
"induction matter" was solved exactly and the 
ID string theory solution was demonstrated. In 
the paper |^] the new critical behaviour of this 
model was claimed, with a log-of-log singularity 



of the string susceptibility. This result contra- 
dicts in a way the results of [|| , so that the whole 
question of critical behaviour in the model needs 
to be clarified, in our opinion. 

In sect. 2 we will formulate the model in the con- 
tinuum and on the lattice, describe the induction 
idea and discuss the planar graphs picture. 

In section 3, the issue of the additional symme- 
try of the centre of the gauge group is reviewed. 

In sect. 4 a general master field equation (MFE) 
will be formulated and the model will be shown 
to be equivalent to a random surface embedded 
in a tree in the large-N approximation. 

In sect. 5 the quantitative approach to the 
model will be formulated, based on its equiva- 
lence to the 2-matrix model (2MM), with the self- 
consistent effective potential for each of matrices. 
Using the orthogonal polynomials approach one 
can then write down a simple functional (and, in 
principle, solvable) equation. As an application 
the induced gauge model with a Gaussian scalar 
potential for any dimension will be solved in sec- 
tion 6 by this method and the result will coinside 
with the known one. 

Sect. 7 will review existing results and proposals 
on the induced gauge theory. 

2. Formulation of the model and its phys- 
ical interpretation 

In the continuum version the model under con- 
sideration is just a gauged scalar field theory in 
the adjoint representation of the group SU(N) 
with the Lagrangian: 

L=^tr ({d,<^ + i{A„<^)f + V{<i>)) (1) 
do ^ ' 

Here $(2;) is an NxN matrix Hermitian field 
(a scalar field in the adjoint representation of 
SU(N)), and A^t,{x) is a gauge field. 

Note that one does not add the Yang-Mills in- 
teraction explicitly. 

The action for the lattice version of the same 
model looks as: 

- E^=i.2,...i,(*(^) - U^{x)^{x + A^)f/t(x))] 
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The partition function is given by: 
yn^"^'*-/ n idU)su{N) exp~NS 



(3) 



<xy> 



The scalar potential can be an arbitrary poly- 
nomial 



(4) 



or even a more complicated function. 

Note that no pure Wilson plaquette term was 
added in @. 

Formally this model can induce the Yang-Mills 
interaction, if one would believe that the param- 
eters of the scalar potential V{(^) could be ad- 
justed so as to get into the regime of asymptotic 
freedom. Of course This is not guaranteed at all 
within this model. Considering this possibility we 
note that already the simplest one-loop logarith- 
mic contribution induces the Yang-Mills interac- 
tion: 

J DAD'i'expi^S) (X J DAexp{-S^nd[U]) (5) 
where the first terms of the 1/m expansion (m^ = 



look as 



1 

yind 



Und[A) = — trF^, 



const 



tr {d,,F^,,f 



.(6) 



with the induced gauge coupling given by 



N 
96^ 



(7) 



1 

dlnd 

From the last formula, we see that the renormal- 
ized mass must be chosen in such a way that 
it would be much smaller than the original cut-off 
scale A, and much larger than the physical mass 
scale jjL of (would be) glueballs. The latter is given 
by the asymptotic freedom relation: 

™2 4g^2 



In 



(8) 



Comparing these two relations we find the scaling 
law: 
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(9) 



(10) 



One should not take this number too seriously. 
This naive scaling law was obtained without tak- 
ing into account the corrections from the scalar 
field self-interactions and from the effects of hard 
gluons (note that higher derivative terms in the 
induced action (|^) are suppressed by l/m^, and 
not by 1/A, and hence will contribute to the 
renormalization of the scaling exponent). 

The correct calculation of this scaling law, as 
well as the whole issue of the existence of the 
correct physical QCD phase in our theory is so 
far beyond our technical possibilities. May be, 
computer simulations, such as those started in 
p^], can shed some light on it. In the opinion 
of the author, it is unlikely that this phase could 
correspond to any local large-N master field ap- 
proach described here, but the possibility for this 
phase to exist for finite N cannot be completely 
excluded. 

3. Planar graph representation and extra 
Zn symmetry 

Let us try to give some geometrical interpreta- 
tion of our theory in terms of planar graph expan- 
sion, since we will try later to give it the meaning 
of some improved bosonic string, and the planar 
graphs usually play the role of a regularization of 
a world sheet of a string theory. 

Let us integrate first over the scalar field in 
(0). In the large N hmit we will get a standard 
planar Feynman diagram technique, as for scalar 
field theory, with <I>'^, -vertices, but the prop- 
agators G^i^i{x, y, A) will be modified owing to the 
external gauge field: they will be given by sums 
over paths of P-ordered Wilson factors in adjoint 
representation: 



G%ix,y,A) 



J DT(x, y) exp —rrfi Length{T) 



(11) 



The whole planar graph expansion can be repre- 
sented in this way as standard planar graphs for a 
scalar field "embedded" into the "external gauge 
field media" . 

The integration that is left, over the gauge field. 
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would produce some non-local interactions of dif- 
ferent pieces of this discretized world sheet. In 
principle, we could try to integrate over them in 
the lattice version of the same representation, us- 
ing the technique proposed in [^HQ, in order 
to get some final representation of this model in 
terms of random surfaces, but it would give any 
promising calculational ideas. 

Let us note only that on the lattice one can im- 
mediately see how the term like the Wilson pla- 
quette action arises in the strong coupling expan- 
sion over the kinetic term in the induced gauge 
theory. Four terms corresponding to the links 
around a plaquette give after a Gaussian scalar 
field integration: 



|tr U {plaquette)Y 



2« trF^^, 



(12) 



As a consequence of the adjoint representation 
used in our model, we will always get the induced 
action for the gauge field with every matrix ele- 
ment of gauge J7-matrix multiplied by some ma- 
trix element of its conjugate . As a result we 
have an extra gauge symmetry in our model, 
the symmetry with respect to the centre of the 
SU{N) gauge group. Namely, one can rotate ev- 
ery U matrix by the Abelian factor: 



.27r 



(13) 



where k^y is an integer, with the consequnce that 
any Wilson average W{C) in the fundamental 
representation will be non-zero (and equal to 1) 
iff the loop C forms a tree in the x-space (with 
zero minimal area of the surfaces spanned on it): 



WiC) l/A^tr;7(C) >= So,AreaiC) 



(14) 



The same will be true for the adjoint Wilson loop 
in the large N limit, since it factorizes in this limit 
into two fundamental loop averages. 

Of course, this kind of superconfinement has 
nothing in common with the physical confinement 
and expected area law for fundamental loops. 

A way out could be some mechanism of spon- 
taneous breaking of this symmetry. Even though 
this is a gauge symmetry, one may hope that it 
may be broken at least in the large N limit. To 
analyse this possibility we can use the arguments 



of yj. For example, for one plaquette action in 
an adjoint representation we can use the large N 
factorization property: 



Sadj — 

Eplaguettes \U {plaquette)\^ 

-N^oo N(3 < j^tc U {plaquette) > 
J2piaquettes tr U (plaquette) 



(15) 



which gives a non-linear equation for the effective 
coupling (3 



(3 = I3W {plaquette, (3) 

The trivial solution 

/3 = 0, W{plaquette, /?) = 



(16) 



(17) 



corresponds to the non-broken Z^r symmetry, 
where as a non-trivial solution 



^ 0, W {plaquette, (3)^0 



(18) 



corresponds to the broken Zjv symmetry with a 
possibility of the correct physical behaviour. 

It is easy to generalize this picture to the whole 
induced action of our model, which will include 
the sums over all loops, and not only over plaque- 
ttes. 

One should stress that the corresponding large 
N phase transition with this symmetry breaking 
should take place before we approach the contin- 
uous limit of the theory. 

This scenario seems to be unlikely for the large 
N master field approach presented below. It 
seems that this centre of the group symmetry re- 
mains unbroken for infinite N. 

4. Master field equation for the eigenval- 
ues of scalar field 

In order to see that our model is exactly solv- 
able in the large N limit let us choose an oppo- 
site order of integration in the functional integral: 
first we will integrate over gauge fields, and then 
over the (eigenvalues of the) scalar field. 

First we demonstrate the idea in the continuum 
version, but a rigorous treatment will be possible 
only on the lattice. 

We introduce, as usual, the "angular" 
parametrization of scalar field in terms of eigen- 
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values 4> = diag{4>i,4>2, ■■■,<t>N) and eigenfunc- 
tions VtijeSU{N): 

^x) ^ VL+ {x)(t>{x)VL(x) (19) 
Putting it into the Lagrangian (0) we obtain: 



where 



B„ 



n+ A^,n{x) + in+ d^,n{x) 



(20) 



(21) 



The field -B^ is just a gauge rotated A^, there- 
fore, naively speaking, the integral over is 
Gaussian. Integrating over it we would obtain 
a simple eff'ective action for only the eigenvalues, 
ready for the application of the large N master 
field approach. But the situation is slightly more 
complicated: any change of variables: PVl, 
(p — > P^(j)P, where P is a permutation matrix, 
does not change $ = rj+(/)J7. Therefore, integrat- 
ing independently over i?^, we would overcount 
the possible configurations of scalar fields. 

To avoid overcounting we can, say, impose the 
condition (pi < (j)2 < ■■■ < (pN, or we have to sub- 
tract one by one the overcounted configurations. 
The second possibility can be naturally realized 
in the lattice version. 

Returning to the lattice let us note, that af- 
ter a change of variables ( |l9| ) at every vertex of 
the hypercubic lattice, we can again choose the 
gauge in such a way that the angular degrees of 
freedom will be absorbed into the gauge fields: 
^xUxy^t ~^ Uxy, because of the group invari- 
ance of the Haar measure (on every link of the 
lattice separately). Then the integral over every 
link gauge variable can be performed by means 
of the so-called Itzykson-Zuber-Kharish-Chandra 
formula p^ : 

/(</), X)^IDU exp (TVtr (0 - UxU^?) 

^ dct.,cxp(Af(0.-Xj)") (22) 



A(0)A(x) 



where 

A(0)-n( 



(23) 



Taking into account extra A^((/)(.t)) of the Dyson 



measure for the " angular" parametrization of the 
scalar field, we arrive at the partition function: 



Z 



(fe) 



exp(-5'e//) 



(24) 



(25) 



where 

Seff = T,xv^OgI{(f>x,<l>y) + 

Ex {^ogA^^x)-Nj:tiV{^'x 

The above mentioned naive continuum action 
( [2^ ) after integration over B^, will correspond 
only to the diagonal term in the determinant in 
([2^). All other (A^! — 1) terms can be considered 
as a consequent subtraction of the overcounted 
permutations of eigenvalues. It would be inter- 
esting to find some continuous field theoretical 
description of this object as a generalization of 
an ordinary spatial derivative. 

Let us remark that in spite the Van-der-Monde 
determinants in the denominator of (E3), the ac- 



tion (25) is not singular at all with respect to 
the coinciding eigenvalues: the determinant in 
the numerator has zeros cancelling these singu- 
larities, since the whole integral evidently has no 
such singularities. 

Now our effective action (|2^ ) depends only on 
N eigenvalues and we are in position to apply 
the saddle-point method in the large N limit and 
derive the corresponding saddle-point equation. 
Namely, we look for the classical configuration for 
the eigenvalues obeying the stationarity condition 
(MFE): 

dSeff 







(26) 



From the physical point of view it is natural to 
expect the spatially homogeneous vacuum for this 
system, so, identifying 



(pl. — (pi — const{x) 

we can write the effective action as 

S'e// = {volu'me){2D)[{\ogI{(p,(p) + 
^^ogA\<P) 

Now the MFE (||) reads as 



(27) 



(28) 



2D- 



2E' 



= V'[(Pk) (29) 
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or, in terms of density of eigenvalues, 



1 dk{cj)) 
N dcf) 



(30) 



the basic quantity to be calculated in this approx- 
imation: 



2D 



d\ogI{ 



2P / d^lp{^^)- 



1 



But we have not yet closed the system of equa- 
tions, since we still have to find some effective 
approach for the calculation of /(</), (p). The orig- 
inal formula of Itzykson and Zuber ( p2[ ) is valid 
for any N, but it is not very suitable in the limit 
of large N, since we have to deal with the sum of 
N! sign-changing terms. 

An approach to this problem was sketched in 
Q , where the integral equation for the calculation 
of the IZ-integral in the large N limit was derived. 
Later we will discuss this approach. 

In the next section we will propose a new ap- 
proach based on the effective 2MM. 

To conclude this section, let us recall a com- 
ment, made in Q, on the equivalence of our in- 
duced gauge model and a multi-matrix model em- 
bedded into an infinite Bethe tree with the coor- 
dination number 2D, in the large N limit. This 
Q-matrix model obeys the following action: 



iVtr 



(32) 



1=1 



where the first sum runs over the bonds of an 
infinite tree with the coordination number 2D, 
and the second over its vertices. Note that if we 
would gauge the model in the spirit of our induced 
gauge theory, it would not change it since, on 
any tree, the gauge variables can be completely 
absorbed into the angular degrees of freedom of 
matrices. 

On the other hand, integrating by means of 
the IZ-integral over the relative "angles" on ev- 
ery bond, and looking for the homogeneous saddle 
point for the eigenvalues (the infinite Bethe tree 
is certainly translational invariant) , we obtain the 
same effective action ( p8| ) and the same MFE ( p9| ) 
or (|3|) as for the D-dimensional induced gauge 
model. 



Hence, in the large N limit the models coin- 
cide. Of course this coincidence is very worrysome 
from the point of view of a possible equivalence 
with QCD. As we see the dependence on the space 
structure almost disappears from the problem in 
the main order in 1 /N. Its dimensionality D enters 
the MFE only through the number of the nearest 
neighbours on the lattice. So, if we would formu- 
late our theory on the 2-dimensional triangular 
lattice, we would obtain an effective dimension- 
ality equal to 3, and for the hexagonal lattice it 
would be 3/2. Of course, in the next 1/N correc- 
tion the dimension of space will enter in a more 
serious way, but we expect even the large N vac- 
uum of QCD to be less trivial. 

On the other hand, none of these features pre- 
vent us from hoping for some meaningful string 
theory for D > 1 following from the induced 
gauge theory. For example, the ID strings ]2^-|2^ 
have an infinite mass spectrum with an appropri- 
ate scaling, and they are a particular case of our 
model. So, we can hope for the same properties 
of it for D>1. 

Another hope (although a weak one) to attack 
QCD starting from our model is the possibility to 
reach the physical QCD behaviour by adjusting 
also N to some critical value which is not sepa- 
rated from A'^ = oo by a phase transition point, 
even though the starting point is not an asymp- 
totically free theory. It resembles the hopes re- 
lated to the strong coupling expansion in Wilson 
lattice QCD which died when confronted with the 
computer simulations. The advantage of the in- 
duced gauge theory is that it gives, at least, most 
probably starting from the zero approximation 
the infinite mass spectrum. 

5. The effective two-matrix model and the 
MFE based on it 

In this section we propose a new approach to 
the induced gauge theory based on its equivalence 
with an effective 2MM. It is completely differ- 
ent from the one presented in our first paper 
which used some 2MM representation for the IZ 
integral. The latter was valid for any 1/n order, 
but the resulting integral equations were difficult 
to analyse. 
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This time we will work directly in the large N 
limit and look from the very beginning for the 
homogeneous MF. In this case, as we know, our 
effective action is given by (|8|). Note that the 
overall factor {volume)2D does not influence the 
MFE (p9|). This fact can be used to reduce the 
problem to an equivalent 2MM with some modi- 
fied effective potential. 

The 2MM action is a particular case of the ac- 
tion (^2|) for two matrices. In terms of two sys- 
tems of eigenvalues {<^fe} and {ipk}, the partition 
function reads: 

Z2MM = /d^0d^Vexp[-iV^f^,(-l/(0fc)+ , . 
y(Vfc))]A2(0)A2(V>)/(0,V) ^ ' 

According to the Z2 symmetry of this model, 
there should exist a symmetric saddle point so- 
lution 0* = ?/)* = Xk ■ 

Comparing the effective action ( ^3| ) with the 
effective action (|8|) at a homogeneous saddle- 
point, one finds that the induced gauge theory 
is equivalent to the 2MM with (the derivative of) 
the effective potential for the latter, which reads 



1 



2D 



2D -I 
DN 



1 



Xk 



-(34) 



It is useful to introduce the density of states, the 
basic quantity to be calculated in the large N 
limit: 



p{x) 



1 dk{x) 



(35) 



N dx 

and the analytical function corresponding to the 
loop amplitude in 2MM: 

1 /■ dypiy) 



W{x) =< tr- 



>= 



(36) 



X ~ ^ J X — y 
The effective potential then reads: 

KV/(x) = ^U'ix) + ^^ReW{x) (37) 

Hence to find W(x) we have to solve a self- 
consistent problem: solving the 2MM with the 
effective potential depending on W(x) we find a 
self-consistency equation on W(x). This trick al- 
lows us to avoid the problem of calculation of 
the IZ integral, since we will now use a well- 
elaborated orthogonal polynomial formalism for 
the 2MM. 



Let us review this formalism in the case of an 
arbitrary potential, using an elegant approach of 
the paper [p^ . 

We start from the partition function of 2MM 
in the form 

Z2MM = 

J d^xd^y/^{x)^{y) (38) 

exp(iVE,(^.yi-n^O-i^(y.))) 

One introduces, according to ||l^, the orthogonal 
polynomials: 

P„(a;) =a;" + a„_ia;"~i + ... + ao (39) 
obeying the orthogonality condition: 

< m\n>= J dxdy 
e^^(^^{xy-V{x)^V{y))) (40) 

Pmix^Pniy) — ^n^mn 

By means of the Jacobi relations 

x\n >= J2i>~i x„,n+i\n - I Xn{z)\n > ^^^^ 
■§^\n >= J2i>iPn.7i+i\n - I >=pn{z)\n > 

where z — exp[— one introduces two opera- 
tor functions which obey the Heisenberg commu- 
tation relation: 



[p„{z),x„{z)] 



1 



(42) 



These relation, can be effectively used in matrix 
models of 2D gravity as a KdV-type approach 
©111. 

In the large N limit we introduce the rescaled 
variables 

= (Art) 

and the KdV-type relation takes the form of the 
Poisson brackets: 



{p{X,uj),x{X,uj)} = 1 



(44) 



where p(A, w), a;(A, w) are already ordinary func- 
tions. We could have started from these KdV- 
type relations, but we will follow the method of 

®- . . 

Considering the matrix element of < m\-^\n > 
and integrating by parts inside it, we obtain the 
equation: 



p{z) = V'{x{z)) - X 



(45) 
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We dropped here the exphcit dependence on the 
cosmological constant A and used the transposi- 
tion operation for x{z) 

x'^{z)=x{^) (46) 



7^' 



= /(A) 



where /„ = 

Eq.(p5|) is in principle sufficient to define both 
X and p. Since we know that p contains only 
positive powers of z, we can compare the coef- 
ficients of non-positive powers. This gives suf- 
ficient number of relations to define everything. 
Particulary simple all this looks for the polyno- 
mial potentials V{x), owing to the finite sums in 
(pi]). One can, say, easily recover all the results 
for the Ising model on random graphs |l^,|T^,^ . 
In a forthcoming paper pT|, we will show how 
effective equations like (H) and (45) are for the 
complete analysis of the critical regimes of 2MM. 

We can now use eq.(p5|) for our D-dimcnsional 
induced gauge theory, where we have now to sub- 
stitute V{x) by Veffix) from iM). We get 



p{z) 



(47) 



^U'{x{z)) + i^^ReW{x{z)) - x{j-J 

To close the equation one has to calculate W{x) 
in terms of x{z). One can obtain this relation 
from an orthogonal polynomial representation of 
the eq. (|3^). ft reads: 

dW{X,x) dlogz{X,x) 

-^T- - ^^^^ 

where z{x) is the function inverse to x{z), to be 
found from (^7|). 

In the next section we will demonstrate how 
this equation works for some particular cases. 

6. Particular examples 

Let us demonstrate how the approach worked 
out in a previous section works for some partic- 
ular examples: for usual 2MM, for _D = (pure 
2d-gravity) and for arbitrary D and quadratic po- 
tential U{x). 

a.Two-Matrix Model: D = 1/2 

In this case Veff{x) = U{x), and we return to 
the original 2MM equation (^). A rather com- 
plete analysis of it will be done in |^ , so we do 
not continue on this subject here. 



b.Pure Gravity: D — 

Now we can retain in (^) only those terms 
thath are singular in l/D, which immediately 
gives the well-known equation for the 1-Matrix 
Model: 



2ReW{x) = U'{a 



(49) 



with all familiar consequences for the pure 2D- 
gravity following from it. 

c. Quadratic Bare Potential for any D 

In this case we have: 



U' {x) = rn^x 



(50) 



It is natural to expect the semi-circle law for the 
distribution of the eigenvalues here, so we will try 
the ansatz: 



W{x) = ^{x- Vx^~4R 
2x1 V 



(51) 



where R{X) has to be found. According to (pq), 
this corresponds to 



x{z) — — h Rz 
z 

and 

ReW{x) = 



(52) 



(53) 



If we insert all these expressions into ( ^ ) we will 
find that p(A, z) should be linear in z and in the 
normalization (M3) it looks like 



p{X, z) — Xz 



(54) 



Comparing the coefficients of ^ and z in (^7|), we 
find: 



1 _ ni'^iD - 1) ± D^m" - 2(2i:> - 1) 



R 



4:{2D- 1) 



(55) 



For > 1 and > m* ^ 2(2D - 1) we choose, 
of corse, the positiv root. This result coincides 
with that one obtained in |^ by Migdal's method 
(see also 0). 

This solution is perfectly valid only in the case 
(which can be called a strong coupling regime) of 
real R or for the effective mass M: 



M'* = - 2{2D - 1) > 



(56) 
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It is clear that the point M = corresponds to the 
critical behaviour of Bethe tree, when its volume 
becomes infinite. It is hard to beleave that some 
interesting (weak coupling) regime exists beyond 
this point. 

d) Upside-Down Quadratic Potential for D=l 
In this case, as we know (see ||2^] for review) 
the critical regime corresponds to "upside-down" 
oscillator potential, so we probably have to take 
the negative root in ( ^ ) and change x into ix in 
(|l|), which would give the right distribution of 
the eigenvalues for the ID-matrix model in the 
critical regime describing the ID-bosonic string 
theory plpil. 



The result can be nicely presented as follows: 



7. Recent developments 

Here we will discuss some alternative ap- 
proaches to the induced gauge theory. 
a.Migdal's Approach 

It is based on an integral equation found in 
defining the large N limit of the Itzykson-Zubcr 
integral. It allowed some new critical behaviour 
to be found for D > 1 in the model, as well as the 
first attempts in the spectrum calculations to be 
made. 

We will look at this equation from a point of 
view that is a bit different of considering it 
as a loop equation for the pure 2MM. Then one 
can easily restore the whole equation for the D- 
dimencional theory, by substituting the bare po- 
tential with the effective one, according to for- 
mula (H). 

Consider an obvious identity for the resolvent 
in 2MM: 



(57) 



W{x) = 
exp (-iVtr y ($) - mr V(l')) 

Integrating by parts in we get 
W{x) = 

J d^^^d^^^expNtT ($* - y($) - y(*)) (58) 



tr 



1 



W{x) = 
i^log(x-r(y) 



Wiy)) 



(59) 



where the contour integral encircles the cuts of 
W{y) but leaves aside the singularity of the log- 
arithm. 

It is the loop equation for the 2MM, which was 
not known before (see for a similar approach) . 

It is sufficient to substitute the bare potential 
with the effective one (|^ in order to recover the 
full D-dimensional equation of [^j. 

Using this approach it was found also in ||] 
that the singularity of W{x) near the end of the 
cut should be x" with a = 1 + ^ arccos( gj^^ ). 
But the a discrepancy with the case D=l/2 (Ising 
model on random graphs) where a = 4/3, is a bit 
confusing to the author. May be this solution is 
aplicable only to 13 > 1. 

b. Log (log) Critical Behaviour 

Another possible solution was advocated in . 
It was found that the system for D > 1 has a crit- 
ical behaviour of effective ID-matrix model with 
the upside-down harmonic potential. The self- 
consistency equation on the "frequency" leads to 



a singularity of the free energy 



ith re- 



The r.h.s. of (58) was calculated in g in the large 
N limit by means of the Ricmann-Cauchi method. 



log log (A) 

spect to the cosmological constant. This is an 
interesting possibility since it would presumably 
give a well-scaled infinite mass spectrum, as for 
the conventional D=l model. Unfortunately it 
can be well justified only in the limit of large D. 

The discrepancy between the two approaches 
is worrysome. It seems that the model still needs 
deeper understanding before we can view it as a 
model of an extended object in physical dimen- 
sions. 

Some related problems and models can be 
found also in |p9| , p7| . 

Let us note also that a similar model in two di- 
mensions, but with an F^^^ term, was investigated 
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